First-principles study of the ferroelectric Aurivillius phase Bi 2 W0 6 
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In order to better understand the reconstructive ferroelectric-paraelectric transition of Bi2W06, 
which is unusual within the Aurivillius family of compounds, we performed first principles- 
calculations of the dielectric and dynamical properties on two possible high-temperature paraelectic 
structures: the monoclinic phase of A2/m symmetry observed experimentally and the tetragonal 
phase of IA/mmm symmetry, common to most Aurivillius phase components. Both paraelectric 
structures exhibit various unstable modes, which after their condensation bring the system toward 
more stable structures of lower symmetry. The calculations confirm that, starting from the para- 
electric A2/m phase at high temperature, the system must undergo a reconstructive transition to 
reach the P2\ab ferroelectric ground state. 



I. INTRODUCTION 

The Aurivillius phases constitute a family of natu- 
rally occurring complex layered oxide structures in which 
fmorite-like (Bi 2 2 ) +2 units alternate with perovskite- 
like (A m _iB m 03 m+ i) -2 blocks, where A= Ca +2 , Ba +2 , 
Sr+ 2 , Bi+ 2 , etc., B = Fe+ 3 , Ti+ 4 , Nb+ 5 , Ta+ 5 , W+ 6 , etc. 
and to is the number of BOe octahedra in the perovskite- 
like blocks (to = 1 — 8) (Tj. The oxidation states of the 
cations and to are interdependent to always provide the 
right number of electrons to oxygen atoms. This fam- 
ily has numerous representatives, the most studied being 
Bi 2 W0 6 (to = 1), SrBi 2 Ta 2 9 (to = 2) and Bi 4 Ti 3 0i 2 
(m = 3) [HO]. 

The ferroelectric properties of the Aurivillius phases 
have been widely studied since the early works of Sub- 
barao [7] and Smolenski et al. [5] and they continue to 
rise technological interest mainly for their potential use 
in thin film non- volatile memory applications [9_. They 
are also attractive at the fundamental level because of 
the peculiarity of their phase transitions and the com- 
plex interplay of their structural instabilities. 

Bi 2 WC>6 is the simplest member of the Aurivillius fam- 
ily with one oxygen octahedron sandwiched between two 
Bi 2 2 layers. It possesses many interesting physical 
properties such as ferroelectricity with large spontaneous 
polarization and high Curie temperature (950°C), piezo- 
electricity with a potential for high-frequency and high- 
temperature applications as an alternative to BaTi03 
and PbZr 1 _ 2 .Ti :I .03 [TU], high ion conductivity [TT] and 
photocatalytic activity |12j . 

Understanding ferroelectricity and the sequence of 
phase transitions in this compound is still challenging. 
Experimental studies [T3] showed that Bi 2 W06 is po- 
lar at ambient temperature with orthorhombic P2\ab 
(No. 29) structure. Rising the temperature reveals at 



670°C a transition to an intermediate polar phase with 
orthorhombic B2cb (No. 41) structure. This structure is 
maintained until 950°C, then the material displays a tran- 
sition to a paraelectric monoclinic A2/m (No. 12) struc- 
ture. This type of ferroelectric-paraelectric transition to 
an A2/m phase is unexpected and rather unique among 
Aurivillius phases, where usually compounds transform 
at high temperature into a centrosymmetric tetragonal 
I4/mmm (No. 139) structure [B]. The same kind of be- 
havior appears in Bi 2 MoC>6, a compound isomorphic to 
Bi 2 W06, which exhibits orthorhombic structures at low 
and intermediate temperatures and a monoclinic P2\/c 
(No. 14) structure at high temperature [141 IT5] . Such 
ferroelectric to paraelectric transitions from orthorhom- 
bic to monoclinic are necessarily reconstructive since the 
space groups of the low temperature phases are not sub- 
groups of the high temperature phases (^42/m or P2i/c), 
but rather of the IA/mmm structure. 

We performed first-principles calculations in order to 
better understand these intriguing observations. First- 
principles studies of Bi 2 WOg have been previously re- 
ported in the literature. Machado et al. [16 performed 
frozen-phonon calculations of the polar modes at the 
Brillouin-zone center in the hypothetical I4/mmm tetrag- 
onal paraelectric phase in order to clarify the microscopic 
origin of the ferroelectric instability. More recently, 
Mohn and Stolen |T7] highlighted the modification of the 
oxygen arrangement around Bi atoms in the Bi 2 2 lay- 
ers at the reconstructive ferroelectric-paraelectric tran- 
sition. In the present work, we complete these stud- 
ies by reporting a systematic investigation of the struc- 
tural, dielectric and dynamical properties of both po- 
tential high-temperature A2/m and I4/mmm paraelec- 
tric phases. We investigate the energy landscape around 
these two phases. The structure and internal energy of 
various other intermediate phases, reached by the con- 
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FIG. 1. The monoclinic A2/m paraelectric phase of F^WOe 
showing edge-sharing octahedra in the ac plane (a) and 
corner-sharing octahedra in the be plane (b) . 



densation of individual or combined lattice instabilities at 
the zone center and zone boundaries of the Brillouin zone, 
are calculated. A diagram comparing the internal energy 
of various possible phases is then established, which con- 
firms that starting from high-temperature paraelectric 
A2/m phase, the system must undergo on cooling a re- 
constructive phase transition to reach the P2\ab ground- 
state structure. The spontaneous polarizations of differ- 
ent polar phases are also reported. 



II. TECHNICALITIES 

The calculations were performed within density func- 
tional theory |181 ITT)] . using a plane waves pseudopo- 
tentials method thanks to the ABINIT package [20H22] . 
The exchange-correlation energy functional was evalu- 
ated within the Local Density Approximation (LDA) as 
parameterized by Perdew and Wang [23^ . Extended norm 
conserving pseudopotentials from Teter [53] were used, 
with Bi(5d 6s 6p), W(5s 5p 5d 6s) and 0(2s 2p) lev- 



TABLE I. Calculated cell parameters and internal atomic 
positions in the paraelectric A2/m phase of E^WOe . Exper- 
imental values from Refll3p are reported between parentheses 
for comparison. 
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els treated as valence states. The wave functions were 
expanded up to a kinetic energy cutoff of 50 Hartrees. 
Integrals over the Brillouin zone were approximated by 
sums on a 6x6x1 Monkhorst-Pack mesh of special k- 
points [25 . The structural optimization was done using 
the Broyden-Fletcher-Goldfarb-Shanno minimization al- 
gorithm (BFGS) [26] - We calculated the ab initio forces 
on the ions and relaxed the position of each individual 
atom until the absolute values of the forces were con- 
verged to less than 10 -5 Ha/Bohr. Phonons, Born effec- 
tive charges and dielectric tensors were calculated using 
Density Functional Perturbation Theory [57] and sponta- 
neous polarization using the Berry phase formalism |29| . 



III. HIGH TEMPERATURE A2/m 
MONOCLINIC STRUCTURE 

We begin our study by fully characterizing the ex- 
perimentally observed paraelectric A2/m phase. In this 
structure for the first time, phonon instabilities are calcu- 
lated together with dielectric properties and Born effec- 
tive charges. Further, we will move on with the character- 
ization of the hypothetical paraelectric 14/mmm phase. 



A. Structure 

The monoclinic A2/m paraelectric phase of Bi2WC>6 
(conventional C2/m) contains 2 formula units per prim- 
itive unit cell (18 atoms). The structure, illustrated in 
Fig.l, consists of alternative Bi2C>2 fluorite-like and WO4 
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TABLE II. Non-zero elements of the calculated Born effective 
charge tensors in cartesian coordinates (see text) for Bi2WOe 
in the monoclinic A2/m phase. The nominal charge of each 
atom (Z nom ) is reported for comparison 
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layers. Within this structure the oxygen octahedra ex- 
hibit a peculiar edge-sharing along a-axis. Along 6-axis, 
the oxygen octahedra are corner-sharing. As conven- 
tional for Aurivillius phases, the stacking axis normal to 
the layers is defined as the c-axis of the A2/m structure. 
The calculated lattice parameters and internal positions 
are shown in Table I. Except for the systematic LDA 
underestimate of lattice parameters, the results are in 
reasonable agreement with experiment. We notice that 
the lattice constants of Mohn and Stolen [T7] within the 
generalized gradient approximation (GGA) are in closer 
agreement with the experiment. This might however be 
partly fortuitous: at the opposite to LDA, GGA is well 
know to overestimate lattice constants and this typical 
error (observed in the lower-temperature B2cb phase, as 
discussed in a next Section) might be partly compensated 
here by the thermal expansion, since the experimental re- 
sults to which the zero-temperature calculations are com- 
pared have been done at 965° C. 



B. Dielectric properties 

The study of the Born effective charges (Z*) and opti- 
cal dielectric tensor (Eoo) is important for the identifica- 
tion of the long-range dipolar contribution to the lattice 
dynamic of a polar insulator. Our results for Z* of the 
non-equivalent atoms are presented in Table II. For clar- 
ity, we note that A2 / m has its a and b axis aligned along 
the cartesian x and y directions, respectively, and c-axis 
makes an angle of (j3 - ?) with z direction. Within this 
choice of cartesian axis, the Born effective charge tensor 
has the form: 
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The Born effective charges for this material are much 
larger than the nominal ionic ones and so follow the ten- 
dency observed for the class of perovskite ABO3 com- 



pounds [3U|. Giant anomalous Z* are observed on W 
(+10. he) and O'i (-7.76e) belonging to the W-O'i chains 
aligned along 6-axis. These values are similar to those 
reported for WO3 [^J and can be assigned to W 5d - O 
2p dynamical changes of hybridization. Anomalous parts 
of the charge on W and on its surrounding 0"i, O^, and 
0"2 oxygen atoms, are strongly reduced along x and z 
directions where W-0 chains are broken. Anomalous Z* 
are also observed in the Bi202 layer with a large contri- 
bution observed on Bi and its surrounding oxygen atoms. 
One can see a significant difference between Bii and Bi2 
(in x direction) despite their apparently similar environ- 
ment in the Bi202 layer. The Born effective charge is 
known to be a sensitive tool to probe covalency effects 
[30] and our results therefore support the recent work of 
Mohn et al. [T7], who suggest that the local environ- 
ment around Bi is asymmetric, with Bii and Bi2 forming 
different number of bonds with the apex oxygen of the 
perovskite layer. 

Finally, the optical dielectric constant is computed: 



/ 6.49 0.08 

7.25 
v 0.08 6.37 



We have no experimental values to compare with and 
our prediction has to be considered as an upper bound, 
because of the typical LDA overestimate of . 



C. Phonon instabilities 

Phonon calculations were performed to study the 
structural stability of the A2/m monoclinic configura- 
tion. The zone-center phonons were computed in the 
conventional cell of 36 atoms, with the primitive cell dou- 
bled along the 6-axis. Due to the folding of the Brillouin 
zone, this gives access to the phonons at T and V. The 
full list of r frequencies is reported in Appendix for com- 
pleteness. The system reveals four instable modes : two 
at r with frequencies at 108z cm -1 and 63i cm -1 , cor- 
responding to the irreducible representations labels [32 
and T%, respectively, and two others at the zone- 
boundary point V=(0, 1/2, 0), with frequencies at 89* 
and 65i cm -1 , corresponding to and , respec- 
tively. The strongest unstable mode is polar along 
6-axis. From Table III, where real-space eigendisplace- 
ments of these unstable modes are given, we can see that 
all of them restrict their atomic motions along 6-axis di- 
rection only and that contrary to what is usually ob- 
served in other paraelectric Aurivillius phases, the A2/m 
monoclinic configuration is not unstable against oxygen 
octahedra rotational modes, which are prevented by the 
edge-sharing spatial organization. 

The unstable modes were then individually condensed 
into the A2/m phase and each of the lower symmetry 
phases reached under these condensations is fully relaxed. 
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TABLE III. Real-space eigendisplacements rj of the unstable phonon modes of the A2/m structure. The eigendisplacements are 
normalized within the 36 atoms cell according to < rj\M\rj >= 1, with the mass matrix M in atomic mass unit |28| . Internal 
coordinates of the atoms of the A2/m phase are those presented in Table I. Modes frequencies between brackets are in cm -1 . 
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TABLE IV. Cell parameters (lengths in A and angles in de- 
grees) and difference of energies AE m with respect to the 
paraelectric monoclinic A/ 2m phase (in meV/formula unit) 
after full atomic relaxation of the different phases reached by 
the condensation of individual unstable phonon modes in the 
A/2m phase. The label and frequency u) ( in cm -1 ) of the 
mode that was initially condensed, together with the symme- 
try of the related phase are also mentioned. 

Mode lo Symmetry Cell parameters AE m 
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From the results given in Table IV, we can see that the 
condensation of the polar unstable mode into A2/m 
leads to a ferroelectric phase with monoclinic A2 (No. 
5) structure, while the condensation of Vjf , V^ and 
lead to monoclinic PI (No. 2) structure. After the full 
relaxation, these phases exhibit nearly identical cell pa- 
rameters and total energy values close to each others. 
The lowest energy configuration corresponds to the mon- 
oclinic PI phase arising from the condensation of Vf 
with a decrease of energy, with respect to A2/m phase, 
of 9.25 meV/formula unit. Further phonon calculations 
in this phase do not reveal any unstable mode, identifying 
it as a potential ground state. 



D. Spontaneous polarization 

As highlighted in the previous section, starting from 
the A2/m phase, a ferroelectric A2 structure can be 
reached from the condensation of an unstable mode, 
giving rise to a polarization along the b direction. The 
amplitude of this polarization can be computed using the 
Berry phase formalism but requires to be cautious. 

Indeed, the modern theory of the polarization (see for 
instance Ref. [29]) teaches us that the spontaneous po- 
larization P s of the ferroelectric A2 structure is related 
to the formal polarizations of the ferroelectric (Pa2) and 
paraelectric (P4.2/771) phases by 



Ps — PA2 - PA1/r, 



n\P< 



Q\ 



(1) 



where n is an integer and Pq = 2 e R/i~2 is the so-called 
polarization quantum, with R the shortest lattice vec- 
tor along the polarization direction and Q the unit-cell 
volume. The factor of 2 accounts for the double band 
occupancy. 

From Eq. [I] the determination of P s from the only 
knowledge of the formal polarizations Pai and PA2/m is 
only possible modulo Pq. This quantum indeterminacy 
arises from the fact that the formal polarizations of the 
Pai and PAijm phases are themselves only defined mod- 
ulo Pq so that the value of P s depends on the specific 
choice of a polarization branch (i.e. choice of how we con- 
nect a given PA2/m point on the left to a given Pai point 
on the right in Figure 2). Of course, P s is a well-defined 
quantity but its unambiguous determination additionally 
requires to determine n through the proper integration 
of the polarization current along an adiabatic insulating 
path between the paraelectric and ferroelectric phases. 

When the polarization quantum is large compared to 
the spontaneous polarization, like in conventional ABO3 
compounds, the uncertainty can be easily fixed. But 
in materials like Aurivillius compounds that have large 
unit cells, the quantum is small and potentially of the 
same order of magnitude as the spontaneous polariza- 



tion (in the present case Pq=25.90 /iC/cm 2 ) and get- 
ting rid of the quantum indeterminacy requires to follow 
the evolution of the polarization along the atomic distor- 
tion from the paraelectric to the ferroelectric phase. In 
practice, this can be achieved by computing the formal 
polarizations of intermediate structures as illustrated in 
Figure 2. The choice of the appropriate branch (solid 
lines in the graph) becomes obvious when the change 
of polarization from one intermediate structure to the 
next one is small compared to the polarization quan- 
tum. From this, we deduce that the spontaneous po- 
larization of the A2 phase is P s — 20.02 /zC/cm 2 . 
We notice that this value differs significantly from the 
one that would have been naively deduced from the dif- 
ference between the values of the formal polarizations 
in the interval [— Pq/2, Pq/2], as typically delivered in 
the output of conventional first-principles softwares (i.e. 
Pa2= -5.88^C/cm 2 and PA2/m=0^C/cm 2 yielding P s = 
-5.88^C/cm 2 ). 

Alternatively, it is also possible to estimate the polar- 
ization from the knowledge of the Born effective charges 
and atomic displacements using the expression: 

Ps > a = o Z taP ST *,P> ( 2 ) 

where 5t K: p is the displacement of atom k along direc- 
tion j3 from the paraelectric to the ferroelectric phase and 
Z* a g the Born effective charge tensor of atom k. This 
provides P s — 20.73 /iC/cm 2 in close agreement with the 
Berry phase estimate and coherent with the quasi-linear 
evolution of the polarization with the structural distor- 
tion seen in Figure 2. We notice that Equation 2 gives 
access to the slope at the origin of the appropriate po- 
larization branch in Figure 2 and provides an alternative 
easy way to identify it. 

IV. HYPOTHETICAL HIGH TEMPERATURE 
IA/mmm TETRAGONAL STRUCTURE 

Starting at high-temperature from the experimentally 
observed paraelectric A2/m structure, one would an- 
ticipate from the previous discussion, a transition to 
a non-ferroelectric PI ground state. Instead of that, 
experiments show that Bi2WC>6 exhibits low tempera- 
ture phases of orthorhombic B2cb and P2\ab symme- 
tries, subgroups of the centrosymmctric IA/mmm parent 
phase, that is the high temperature paraelectric phase of 
all Aurivillius compounds with a number of perovskite 
blocks m>l. We so continue our study by considering 
the hypothetical IA/mmm paraelectric phase. 

A. Structure 

The IA/mmm structure as illustrated in Fig. 3(a), con- 
tains 2 formula units per primitive unit cell (18 atoms) 
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FIG. 2. Change in polarization P along a path of atomic 
distortions from the paraelectric A2/m structure to the fer- 
roelectric A2 structure. The possible values of P at given 
atomic positions differ by multiples of the polarization quan- 
tum, Pq=25.90 /zC/cm 2 . Squares represent formal polariza- 
tions in the interval [— Pq/2, P<5/2](grey area). 



and consists of Bi2C>2 fluorite-like layers alternating with 
perovskite-like blocks of corner-sharing octahedra chains. 
Since IA/mmm is not observed experimentally, an esti- 
mate of the paraelectric unit cell is deduced from the 
orthorhombic lattice parameters of the experimental fer- 
roelectric phase cited in Ref. |16| . using a tetra = b tetra = 
3.86,4 (= ^Jaorthborth/^) and c tetra = 16. 52,4 . In order 
to facilitate the comparison with the experimentally ob- 
served orthorhombic ferroelectric ground state, that con- 
tains 4 formula units per primitive unit cell (36 atoms), 
our calculations were performed in a tetragonal IA/mmm 
unit cell doubled along the [110] directions, where the 
lattice parameters become a = b = 5.45 A, c = 16.52A 
(see Fig. 3 (b)), allowing the paraelectric and ferroelectric 
phases to have the same number of atoms and the same 
orientation. In what follows, a,b and c will be the vectors 
of our doubled unit cell so that a and b directions will re- 
fer to the [110] directions of the conventional centered 
unit cell and c to the direction perpendicular to them. 

The structural parameters and internal atomic posi- 
tions of the IA/mmm hypothetical structure after full 
atomic relaxation are reported in the top part of Ta- 
ble V. Our internal coordinates show good agreement 
with those calculated by Machado et al. in Ref.[TB] who 
worked however at larger lattice constants, correspond- 
ing to experimental values estimated as discussed above. 

We further notice that the internal energy of the re- 
laxed IA/mmm structure is 156 meV per formula unit 
higher than that of the A2/m phase. From now, we will 
report the energies taking the IA/mmm structure as ref- 
erence. A summary of the internal energies of all the 
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FIG. 3. The hypothetical tetragonal I4/mmm paraelectric 
phase of Bi2WOe: side view (a) and top view (b), highlighting 
the corresponding doubled cell in [110] direction. 



phases studied will be provided in Section V. 



B. Dielectric Properties 

Our results for the dynamical effective charges of 
Bi2WC>6 in the IA/mmm phase, for non-equivalent atoms 
are presented in Table VI. In this case, the cartesian x, 
y, and z axis are aligned along the crystalline a, b and c 
directions. 

As expected, there are large anomalous contributions 
to Z*, on W(+lle) and Oi(-4.88e) atoms in the octahe- 
dral afr-plane and on W(+9.32e) and C>2(-4.74e) atoms 
in the apical direction. Within Bi202 layers, the large 
charges also observed on Bi(+5e) attest of the ionico- 
covalent character for bonding in both the fiuorite-like 
and perovskite-like units of Aurivillius structures. Com- 
parison with A2/m shows that Z* on W along W-0 
chains are in the same range of values for both structures, 



TABLE V. Calculated unit cell parameters and internal 
atomic positions in the hypothetical paraelectric I4/mmm 
phase doubled along [110] and the B2cb and P2±ab ferroelec- 
tric phases. Experimental values estimated from ref. |16| for 
the tetragonal paraelectric phase (see text) and taken from 
ref. [13 for orthorhombic ferroelectric phases are reported in 
parentheses for comparison. We are not reporting Wyckoff 
positions for the I4/mmm phase since our unit cell is dou- 
bled. 



Structure 


Atom 


Wyckoff' 


X 


y 


z 


IA/mram 












a=5.3ll 


Bi 




0.0000 


0.5000 


0.1756 


(5.45) 






(0.0000 


0.5000 


0.1701) 




W 




0.0000 


0.0000 


0.0000 


c=16.14A 


Oi 




0.2500 


0.7500 


0.0000 


(16.52) 


0-2 




0.0000 


0.0000 


0.1143 








(0.0000 


0.0000 


0.1152) 




o 3 




0.2500 


0.2500 


0.7500 


B2cb 












a=5.35i 


Bi 


8b 


0.9961 


0.5149 


0.1702 


(5.49) 






(0.9947 


0.5116 


0.1722) 




W 


4a 


0.0000 


0.0000 


0.0000 


b=5.38A 


Oi 


8b 


0.2173 


0.7534 


0.9828 


(5.52) 






(0.2981 


0.7399 


0.9892) 




o 2 


8b 


0.0641 


0.9536 


0.1101 


c=16.12A 






(0.0753 


0.9587 


0.1087) 


(16.54) 


o 3 


8b 


0.2369 


0.2547 


0.7484 








(0.2771 


0.2550 


0.7499) 


P2 1 ab 












a = 5.30A 


Bii 


4a 


0.0109 


0.4843 


0.1698 


(5.45) 






(-0.0126 


0.5113 


0.1730) 




Bi 2 


4a 


0.0036 


0.5074 


-0.1698 


b=5.32y4 






(-0.0113 


0.4762 


-0.1719) 


(5.48) 


W 


4a 


0.0000 


0.0000 


0.0000 




O'i 


4a 


0.2318 


0.6896 


0.0108 


c=16.17i 






(0.2679 


0.6933 


-0.0147) 


(16.47) 


0"i 


4a 


0.3321 


0.1984 


-0.0109 








(0.3342 


0.2220 


0.0163) 




0' 2 


4a 


-0.0684 


-0.0368 


0.1109 








(0.0854 


0.0449 


0.1072) 




0" 2 


4a 


0.4428 


0.4563 


-0.1108 








(0.5703 


0.5526 


-0.1086) 




0' 3 


4a 


0.2676 


0.2413 


0.7515 








(0.2728 


0.2508 


0.7489) 




0" 3 


4a 


0.2752 


0.2465 


0.2516 








(0.2740 


0.2326 


0.2514) 



but O'i (-7.76e) in A2/m is much larger than Oi(-4.88e) 
in IA/mmm. 
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TABLE VI. Non-zero elements of the calculated Born effective 
charges tensors in cartesian coordinates for Bi2WC>6 in the 
hypothetical tetragonal I4/mmm symmetry. 



Atom 


7* 

^xx 


7* 
^yy 




Znom 


Bi 


5.00 


5.00 


4.29 


+3 


W 


11.09 


11.09 


9.32 


+6 


Oi 


-4.88 


-4.88 


-1.46 


-2 


2 


-2.7 


-2.7 


-4.74 


-2 


0, 


-2.95 


-2.95 


-2.74 


-2 



Optical dielectric tensor is given by: 

/ 7.60 

£00 = 

V 



7.60 
6.84 



These values are comparable to those reported for the 
A2/m phase. 



C. Phonons Instabilities 

Starting from fully relaxed tetragonal prototype struc- 
ture, we computed the zone-center phonon modes of the 
doubled 36 atoms unit cell that, due to the folding of 
the Brillouin zone, combine in fact phonons at the T 
and X points of the Brillouin zone of the primitive unit 
cell. We identified five unstable phonon modes, twofold 
degenerated along a and 6 directions. The mode with 
the largest instability lies at the Brillouin zone center 
r with a frequency at 198z cm -1 and corresponding to 
the irreducible representation . The four other unsta- 
ble modes are at the zone boundary point X=(~, ~, 0) 
|33| with frequencies at 183« ,135i, 104i and 98i cm -1 , 
and corresponding to the irreducible representations, 
■ X2 , X£ and X£ , respectively. The latter four modes 
are antiferrodistortive modes since they are at the zone 
boundary X; their condensation in the tetragonal struc- 
ture will necessarily lead to the unit cell doubling along 
the [110] direction. The complete list of T phonons and 
a comparison with the results of Machado et al. [16 is 
provided in Appendix. 

In what follows, we will describe the atomic motions 
associated to the different unstable modes as deduced 
from the inspection of the real-space eigendisplacemcnts 
are reported in Table VII. 

The strongest unstable mode (lj = 198i cm -1 ), de- 
generated into r^~(a) and r^~(b), corresponds to the polar 
motion of W and Bi atoms against the oxygens along the 
a-axis for (a) or along 6-axis for (b) , together with 
a small additional shift of Oi octahedral oxygen atoms 
in the a6-plane. 

We notice that, in our calculations, only one polar 
mode is unstable, while the presence of a second un- 
stable polar mode, related to a nearly rigid motion of 



Bi2C>2 layer relative to the perovskite-like block (rigid- 
layer mode), was previously reported in Ref. 1161 Such 
a polar mode exists in our phonon calculations but ap- 
pears at a slightly stable frequency of 29 cm -1 (labeled 
as 17 [29] in Table VII.). 

This difference probably arises from the slightly differ- 
ent cell parameters used in both calculations. 

The second unstable mode X^ (u = 183i cm -1 ), de- 
generated into X2 (a) and X^ (b) , corresponds to an anti- 
polar motion with displacements of the atoms very simi- 
lar to the mode, except for O3 atoms. 

The third unstable mode X^ (oj = 135z cm -1 ), degen- 
erated into X% (a) and A^~(b), consists in a rotation of 
oxygen octahedra about the c-axis. In this mode, W and 
O2 atoms have no displacements while Oi atoms shift in 
the a6-plane. The difference between A"^~(a) and (b) 
lies in a slight motion of O3 along a-axis in (a) or 



The fourth unstable mode X^ (ui = 104? cm 1 ), de- 



along 6-axis for X£ (b). 

The fourth unstable 
generated into A^~(a) and X^(b), consists essentially in 
a tilting of oxygen octahedra around a-axis for AT+(a) 
and around 6-axis for X^(b). In this mode there is no 
displacement for W atoms while Oi atoms move along 
the c-axis. O2 atoms move along the 6-axis for X^(a) 
or along a-axis for X^(b). In the Bi2C>2 layer only O3 
atoms are moving along c-axis. 

The fifth unstable mode X± (ui = 98i cm -1 ), degener- 
ated into X^(a) and X^(b) is very similar to X^ except 
for the motionless O3 atoms and the displacement of Bi 
atoms in the Bi202 layer. 



D. Intermediate and low temperature phases 

The condensation of these different unstable modes 
within the paraelectric I4/mmm phase will lower the 
symmetry and decrease the energy. We have seen pre- 
viously that the internal energy of the IA/mmm phase is 
significantly higher than that of the A2/m phase but its 
phonon instabilities are also about twice larger. This let 
us anticipate that the condensation of the latter could 
produce a decrease of energy one order of magnitude 
larger than those reported in the A2/m phase [34 and so 
yield phases with an internal energy eventually lower that 
that of the PI phase discussed in the previous section. 
The symmetry resulting from the condensation of indi- 
vidual unstable modes in the IA/mmm phase, together 
with the decrease in energy with respect to IA/mmm 
after full atomic relaxation within this given symmetry, 
are reported in the top part of Table VIII. Although the 
different unstable modes are twofold degenerated, we re- 
stricted ourselves to the condensation of one of the two 
degenerated modes since other combinations do not ap- 
pear in the low-symmetry phase experimentally observed 
that are discussed later. 

The condensation of the polar T^"(a) mode in the 
tetragonal paraelectric IA/mmm phase distorts the 
structure into a ferroelectric Fmm2 (No. 42) phase lead- 
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TABLE VII. Real-space eigendisplacements r\ of unstable and some relevant stable phonon modes of the hypothetical 14/mmm 
structure. The eigendisplacements are normalized for the 36 atoms cell according to < r/\AI\rj >= 1, with the mass matrix M in 
atomic mass unit |28| . The internal coordinates of the different atoms are those of Table V. The phonon frequencies, reported 
in brackets, are in cm -1 . 



Atom q 







r 5 "(a) 






X- (a) 






X+(a) 








[198t] 






[183i] 






[135i] 






X 


y 


z 


X 


y 


z 


X 


y 


z 


Bi 


0.0020 


0.0000 


0.0000 


0.0040 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


W 


0.0170 


0.0000 


0.0000 


0.0190 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


Oi 


-0.0400 


-0.0020 


0.0000 


-0.0430 


-0.0050 


0.0000 


-0.0620 


-0.0620 


0.0000 


o 2 


-0.0660 


0.0000 


0.0000 


-0.0610 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


o 3 


-0.0150 


0.0000 


0.0000 


0.0000 


-0.0010 


0.0000 


0.0002 


0.0000 


0.0000 



X+(a) 
[104i] 





X 


y 


z 


X 


Bi 


0.0000 


0.0000 


0.0000 


0.0000 


W 


0.0000 


0.0000 


0.0000 


0.0000 


Oi 


0.0000 


0.0000 


0.0470 


0.0000 


o 2 


0.0000 


0.0640 


0.0000 


0.0000 


o 3 


0.0000 


0.0000 


-0.0060 


0.0000 



*4 + M 






17 




|98i] 






[29] 




y 


z 


X 


y 


z 


0.0060 


0.0000 


0.0150 


0.0000 


0.0000 


0.0000 


0.0000 


-0.0270 


0.0000 


0.0000 


0.0000 


-0.0510 


-0.0230 


-0.0010 


0.0000 


0.0690 


0.0000 


-0.0200 


0.0000 


0.0000 


0.0000 


0.0000 


0.0010 


0.0000 


0.0000 



X 2 X t r 5 







[80] 






[107] 






[137] 






X 


y 


z 


X 


y 


z 


X 


y 


z 


Bi 


0.0220 


0.0000 


0.0000 


0.0000 


0.0220 


0.0000 


0.0010 


0.0000 


0.0000 


W 


0.0050 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


-0.0040 


0.0000 


0.0000 


Oi 


-0.0300 


0.0110 


0.0000 


0.0000 


0.0000 


0.0270 


0.0630 


0.0420 


0.0000 


o 2 


0.0070 


0.0000 


0.0000 


0.0000 


0.0300 


0.0000 


-0.0420 


0.0000 


0.0000 


o 3 


0.0000 


-0.0050 


0.0000 


0.0000 


0.0000 


-0.0080 


-0.0130 


0.0000 


0.0000 



TABLE VIII. Differences of energies with respect to the hypothetical paraelectric tetragonal M/mmm phase, AE t (in 
meV/formula unit), after full structural relaxations of the different phases arising from the freezing of one or more unsta- 
ble modes. The space group of each phase is given in non conventional setting with c set as the long axis. The labels and 
frequencies u (in cm -1 ) of relevant phonon modes are also provided. 



Mode label 




Symmetry 


AE t 


Phases arising from single mode condensation 








17(a) 


198i 


Fmm2 


-114.75 


X 2 -(a) 


183i 


Amma 


-92.52 


X+(a) 


135i 


Abam 


-99.63 


X+( & ) 


104z 


Bmab 


-71.01 


X+(a) 


98i 


Amaa 


-50.04 


Experimentally observed phases 








r 5 -(a)+ X 3 +(a) 




B2cb 


-187.02 


T-(a) + X 3 +(a) +X+(b) 




P2 1 ab 


-216.27 


Hypothetical phases 








r 5 "(a) + X+(b) 




A2\ma 


-206.11 


X 3+ (a) + A 2+ (b) 




P2 1 /a 


-105.43 


r 5 ~(a) +X 3+ (a) +X 2+ (a) 




Pa 


-213.05 


r 5 -(a) +X 3+ (b) +X 2+ (b) 




P2i 


-204.21 
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ing, after full structural relaxation, to cell parameters 
a=5.37i, b=5.36i and c=16.00l. 

The condensation of the X^ (a) mode leads to Amma 
(No. 63) symmetry (conventional CmCm) and relaxed 
cell parameters a=5.36A b=5.34A and c=16.00A 

The condensation of X% (a) mode leads to a or- 
thorhombic structure Abam (No. 64) (conventional 
Cmca) with quasi-tetragonal cell parameters a sa b = 
5.23A and c = 16.361. 

The condensation of X^ (a) mode leads to an or- 
thorhombic structure Bmab (No. 64) and relaxed cell 
parameters a=5.28A, b=5.32A and c=16.26A 

Finally, the condensation of (a) leads to Amaa (No. 
66) symmetry (conventional Cccm) and relaxed cell pa- 
rameters a=5.28A b=5.3L4 and c=16.24A 

We observed that, amongst these phases, the largest 
decrease of energy arises from the condensation of the 
largest instability. The related Fmm2 phase does 
however still exhibit phonon instabilities. 

The ferroelectric intermediate B2cb phase, observed 
experimentally, results from the combined condensation 
of two unstable modes T^"(a) and A^ (a) in the paraelec- 
tric IA/mmm phase, while the addition of a third mode 
X^~(b) allows the ferroelectric P2\ab ground state to be 
reached. We notice that a proper (a) or (b) orientation of 
the unstable modes is necessary to match the experimen- 
tal structure |16| . with a conventional a-axis direction for 
the polarization. 

The relaxed structural parameters of B2cb and P2\ab 
phases are presented in Table V and show the same good 
agreement with experimental data as the other phases 
previously discussed. For the B2cb phase we notice that, 
while LDA underestimates the lattice constants, in this 
case, the GGA calculation of Mohn and Stolen [T7] over- 
estimates them clearly. The related energy, with respect 
to the IA/mmm phase, are reported in Table VIII. The 
lowest internal energy is associated to the P2iab phase, 
which does moreover not exhibit any further phonon in- 
stability, in agreement with the fact that it corresponds 
to the experimental ground state. 

We notice that there are many more ways to combine 
unstable phonon modes than those producing the B2cb 
and P2\ab phases and we have explored few of them. We 
notice for instance that the hypothetical A2\ma phase 
(No. 26) resulting from the combined condensation of 
r^"(a) and A.^ (b) mode is never observed experimentally, 
although it has a lower internal energy than the observed 
intermediate B2cb phase. The decrease of energy associ- 
ated to the two phases B2cb and A2\ma is close to the 
sum of decreases of energy obtained by each individual 
instability producing them, emphasizing the weak inter- 
actions between the rotations and the polar mode. At 
the opposite, the P2\/a phase (No. 14) arising from the 
condensation of X^ (a) and X^ (b) is much higher in en- 
ergy emphasizing the fact that both types of rotations 
tend to exclude each others but not totally. Their com- 
bination with the polar mode will so produces the P2^ab 
ground state. 



At the end of Table VIII, we also check the sensitivity 
of the energy to the relative orientation of the different 
unstable modes appearing in the P2\ab ground state. We 
see that replacing the A^~ (a) mode by Ag" (b) to produce 
a P2\ phase (No. 4) has a larger effect on the energy, 
than replacing the X^/ih) mode by A.^ (a) to produce 
a Pa phase (No. 7). We notice that both A^(a) and 
X^"(b) modes correspond to oxygen rotations about the 
c-axis and only differ by slight O3 motions along either 
a- or b-axis. 



E. Coupling of lattice modes and energy landscape 

Recently, the coupling of phonon modes in perovskite 
compounds has known a significant renewal of interest. 
For instance, the trilinear coupling between the polariza- 
tion and two rotational modes in layered perovskite, like 
artificial (SrTi03)/(PbTiC>3) superlattices [35] or natu- 
ral Ruddlesden-Popper Ca3Ti2C>7 compounds (35], has 
been identified as a promising road to achieve new and/or 
enhanced functional properties. Aurivillius phases like 
SrBi2Ta20g are known to exhibit a similar type of tri- 
linear coupling |3J and it is interesting to question if the 
same type of phenomena is also inherent to Bi 2 WOq . In 
what follows, we further quantify the contributions of in- 
dividual phonon modes to the different phases previously 
discussed in order to get a better description of the en- 
ergy landscape and significant phonon couplings. 

The set of phonon eigendisplacements rji of the the 
IA/mmm phase (normalized such that < 77 1 1 77 >= 
1)|27 defines a complete basis for the atomic distortions 
from IA/mmm to a given phase, X, so that we can de- 
compose the atomic distortion vector A from IA/mmm 
to A as a sum of contributions from the different modes 

A = J2QiVi (3) 

i 

where the mode amplitude Qi is defined as 

Qt = Aon =< r)i\M\A > . (4) 

The mode amplitude can be further decomposed into the 
distortion amplitude A =< Ai|M|A > and the cosine di- 
rector cti with respect to distortion direction m such that 
^A af = 1. The mode amplitude Qi quantifies which am- 
plitude of rji has been condensed into IA/mmm to reach 
X, while on better describes the relative contribution of 
the different modes entering into a given distortion. 

In the top part of Table IX, we report the most signifi- 
cant mode contributions to the atomic distortions in the 
phases arising from the condensation of individual insta- 
bilities (limiting ourselves to contributions on > 0.01). 
We can see that the Abam and Amaa phases arise from 
the condensation of the related unstable mode only, while 
Fmm2 , Amma and Bmab phases also involve small con- 
tributions of additional stable modes of the same symme- 
try, as induced by phonon-phonon anharmonic couplings. 
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TABLE IX. Contributions 014 (see text) of the different phonon modes i to the atomic distortion from the paraelectric IA/mmm 
phase to the different intermediate and ground state ferroelectric phases of Bi2W06- The distortion amplitude A (Bohr) is 
also reported in each case. Each phonon mode is identified by its symmetry label and frequency (in brackets, in cm -1 ). The 
corresponding real-space eigendisplacements are reported in Table VII. We limited ourselves to contributions a.i > 0.01. 



Phase 


A 


IT 










IT 


X2 


xt 


r 5 " 






[l98i] 


[l83i] 


[l35i] 


[l04i] 


[mi] 


[29] 


[80] 


[107] 


[137] 


Phases arising from single mode condensation 






















Fmm2 


361.0 


0.84 










0.51 






0.11 


Amma 


395.6 




0.79 










0.60 






Abam 


444.7 






0.99 














Bmab 


444.5 








0.98 








0.14 




Amaa 


384.0 










0.99 










Experimentally observed phases 






















B2cb 


601.0 


0.53 






0.79 




0.24 






0.09 


P2 ia b 


582.7 


0.52 




0.59 


0.55 




0.04 




0.1 


0.09 


Hypothetical phases 






















A2\ma 


514.5 


0.58 




0.78 






0.01 






0.11 


P2i/a 


455.0 






0.76 


0.62 








0.11 





As highlighted in the previous section, the ferroelec- 
tric P2\ab ground state results from the condensation of 
three distinct atomic distortions, , X$ and X£ into 
IA/mmm. Due to symmetry considerations, these three 
modes cannot couple linearly. However, the Bi2W06 
structure does allow numerous trilinear coupling terms. 
Thanks to ISOTROPY software [33], we have identi- 
fied the following symmetry-allowed possibilities involv- 
ing the r^" mode and one of the X^ and X^ modes : 

© x+ © A7 
© x+ © xi 

® X+ © ATf 

© x+ © x^ 

This teaches us that the presence of and X% major 
modes in the intermediate B2cb phase, automatically al- 
low the appearance of minor Xj~ and X^ modes. More- 
over the addition of the X^ mode in the P2\ab ground 
state is compatible with the additional appearance of mi- 
nor Xj~ and X^~ modes. If we look in Table X however, 
to the modes that have been condensed to reach the B2cb 
and P2\ab phases, we only see significant contributions 
from r^, X^ and X£ modes (the unstable ones and, 
eventually, higher-frequency modes of the same symme- 
tries). This means that, although various trilinear cou- 
pling terms exist, none of them is strong enough to suc- 
ceed condensing significantly any of the minor modes (i.e. 
ol% is not rigorously zero but smaller than 0.01). 

Since the contributions of all minor modes are neg- 
ligible, we performed an expansion of the energy (in 
meV/formula unit) limiting ourselves to the major unsta- 
ble modes that enter the observed orthorhombic phases 



only (i.e. 1^ (a), X 3 +(a) and X+(b)) : 
AE t (Q r - , Q x + , Q x + ) = 

- 2.50Q 2 - 7.50 x lO" 1 ^ + 

5 3 

- 1.03Q^+ + 1.36 x 1(T B Q*- 

+ 1.97 x 10 -6 Qi+ + 2.65 x 10 -6 Qi+ 
+ 2.10 x lQ- 6 Q 2 r -Q 2 x + 

- 1.79 x 1Q- 9 Q 2 T -Q 2 X + 
+ 1.79 x 10- 6 Q 2 + Q 2 + 

Q represents the mode amplitude as explained in the pre- 
vious section . The values of the polynomial coefficients 
have been obtained by systematic energy calculations af- 
ter the condensation of various combination of the un- 
stable modes. 

The coefficients of Q 2 _ , Q 2 + and Q 2 + are negative. 

1 5 X 3 X 2 

coherently with the phonon instabilities. The biquadratic 
coupling coefficient between rg~(a) and X% (b) is nega- 
tive but its value is extremely small. The biquadratic 
coupling coefficients between X% (a) and both (a) and 
X£ (b) are positive but still rather small, meaning that 
although the appearance of X^ (a) tends to exclude the 
two other mode, this effect is rather small. This is con- 
firmed by the fact that, although a r - decreases due to 
the appearance of the other modes, the mode amplitude 
Q r - remains very similar in Fmm2 1 B2cb and P2\ab 
phases. 

The coupling scenario in Bi 2 WOg is therefore rather 
different from those of other previously studied Aurivil- 
lius phases, namely, Srf^T^Og (777=2) and Bi4Ti30i2 
(to=3). In SrBi 2 Ta 2 0g, the strongest instability is the 
octahedral tilting mode X^ around a/b direction and 
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TABLE X. Spontaneous polarization (in /iC/cm 2 ) of different 
observed or hypothetical polar phases of E^WOe. P S (BP) 
refers to the polarization computed using the Berry Phase 
approach (Eq.l). P 3 (Z*) corresponds to the polarization de- 
duced from the knowledge of the Born effective charges and 
atomic displacements (Eq.2). The latter has been decom- 
posed into the contributions coming from the lowest frequency 
T$ polar modes [57] , identified by their frequency Q (in brack- 
ets, in cm -1 ) and mode polarity S ~}2 Z* rj (in parenthesis, 
in /zC/cm 2 ). 



Phase 


17 


r 5 " 


r s " 


Ps{Z*) 


Ps(BP) 




[l98i] 


[29] 


[137] 








(0.18) 


(0.03) 


(-0.06) 






Fmm2 


55.26 


4.82 


-2.40 


58.59 


53.98 


B2cb 


57.92 


3.77 


-3.30 


58.29 


52.06 


P2 1 ab 


55.09 


0.61 


-3.17 


52.63 


48.16 


A2\ma 


53.32 


-0.13 


-3.42 


49.96 


46.98 



the condensation of this mode alone results in a nonpo- 
lar intermediate phase. This tilting mode and the polar 
mode r^" are mutually exclusive, and the final ferroelec- 
tric phase Ama2i needs the presence of a third stable 
mode X2 to stabilize [3J. In Bi4Ti30i2, as for Bi2W06, 
there is three major unstable modes r^~, X£ and X% '. 
However, there is no occurrence of experimentally de- 
tected intermediate phases, and the ferroelectric JSlal 
monoclinic ground state is reached by the simultaneous 
condensation of the three unstable modes and could be 
denoted as an avalanche phase transition [3]. 

F. Spontaneous polarization 

In both orthorhombic ferroelectric phases, the spon- 
taneous polarization P s expands along a-direction only, 
with no component along c. We report, in Table X, 
the spontaneous polarizations of different observed and 
hypothetical phases as obtained using the Berry phase 
approach (P S (BP)) and using Eq.(2) and the knowl- 
edge of the Born effective charges and atomic distortions 
(P S (Z*)). Again, determination of P S (BP) requires care- 
ful determination of the polarization quantum (see dis- 
cussion in Section III-D). 

In Table X we also report the contributions to the full 
polarization arising from the individual polar modes. 
We see that the main contribution arises from the unsta- 
ble mode, which is also the most polar. We observe 
that the spontaneous polarization is rather constant in 
the different phases, in agreement with the weak cou- 
pling between polar and rotational modes, yielding an 
almost constant amplitude of Q r - as discussed in the 
previous Section. The polar mode of highest frequency 
has a tiny contribution but tends to reduce the full po- 
larization. We notice that P S (Z*) slightly overestimates 
P S (BP) , but since the same distortion is used in both 
approaches this has to be linked to the evolution of the 



Born effective charges along the path of atomic distortion 
and neglected in the computation of P S (Z*). 

Independently, we also decomposed the polarization in 
terms of contributions coming from the individual layers: 
the main contribution (76 % of total polarization) comes 
from the perovskite-like block, due to the displacement of 
W opposite to its surrounding octahedral oxygen atoms 
Oi and O2 while a smaller part (24 % of total polariza- 
tion) comes from the Bi2C>2 layer, due to Bi displace- 
ment opposite to O3. This means that the contribution 
of Bi 2 02 to the polarization is not dominant, contrary to 
the usual belief [B]. 

V. DISCUSSION 

In the previous Sections, we have discussed many dif- 
ferent phases of Bi 2 W06- In Figure 4, we provide a sum- 
mary comparing the internal energies of all the phases 
we have studied. This diagram gives a global view and 
helps clarifying the phase transitions. 

At high temperature Bi2WOg crystallizes in a mon- 
oclinic A2/m structure. Starting from this configu- 
ration, displacive phase transitions are possible to A2 
or PI phases. However, Bi 2 WOg can adopt another 
polymorphic IA/mmm tetragonal form, which displays 
much stronger phonon instabilities, yielding orthorhom- 
bic B2cb and P2iab phases of lower energies. Starting at 
high temperature from the A2/m structure, the system 
has to exhibit a reconstructive phase transition to reach 
the P2\ab ground-state. 

VI. CONCLUSIONS 

We have performed a first-principles study of Bi2WOg 
using Density Functional Theory within the Local Den- 
sity Approximation. 

First, we characterized the high-temperature paraelec- 
tric monoclinic A2/m phase. Phonon frequencies and 
dielectric properties were computed. This monoclinic 
structure, with a configuration mixing edge-sharing and 
corner-sharing (WC>4)~ 2 octahedra, exhibits like in typ- 
ical perovskites anomalous Z* on W and its surround- 
ing oxygens along the corner-sharing chains, yielding po- 
tentially a non-negligible spontaneous polarization along 
this direction. Lattice dynamics study reveals four struc- 
tural insabilities that lead the system to energetically 
quasi-equivalent monoclinic phases. One of these phases 
is ferroelectric with space group A2 and a spontaneous 
polarization of 20 /iC/cm 2 . 

Second, we considered a hypothetical paraelectric 
IA/mmm phase that, although not observed experimen- 
tally, has a symmetry that is the supergroup of the low 
temperature orthorhombic phases, observed experimen- 
tally. We fully characterized the IA/mmm phase, that 
exhibits five degenerated phonon instabilities. Our cal- 
culations agree with the observation that Bi2WOg has a 
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FIG. 4. Summary Diagram of phase transitions in IfeWOe. AE t is in meV/formula unit and represents the difference of energy 
with respect to IA/mmm taken as the zero energy reference. 



P2±ab ground-state structure. This ground state can be 
reached from the condensation of three unstable phonons 
in the paraelectric IA/mmm phase : a polar mode, rj, 
a rotation of oxygen octahedra around c-axis, and 
a tilt of oxygen octahedra around a-axis, X+. The ex- 
perimentally observed B2cb phase appears as an inter- 
mediate phase along this path, arising from the conden- 
sation of the r^" and modes. Our calculations reveal 
a weak coupling between the different modes entering 
the ground-state and the absence of significant distortion 
arising from minor modes allowed by symmetry. As in 
other perovskite ferroelectric oxides, the Z* are anoma- 
lously large, yielding a strong spontaneous polarization 
of 48 /iC/cm 2 in the P2\ab phase. The significantly lower 
internal energy of the Pl^ab phase, compared to the PI 
phase, is in agreement with the fact that starting from 
A2/m phase, the system displays a reconstructive phase 
transition. 
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VIII. APPENDIX 

In this Section we report for completeness the full set 
of r phonons in the A2/m (see Table XI) and IA/mmm 
phases (see Table XII) of Bi2WOg as computed within 
the LDA in the fully relaxed structures. For the A2/m 
phases no previous calculations have been reported and 
no experimental data are available. Our results therefore 
constitute a theoretical prediction that might be useful 
to interpret future experimental measurements. For the 
hypothetical lA/mmm phase, our results are compared 
to those of Machado et al. [TB] also within the LDA but 
at a larger volume. As usual in ferroelectric oxides, the 
use of a larger volume favors the ferroelectric instability. 
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TABLE XI. Frequencies (in cm x ) of the phonon modes at F in the monoclinic A2/m phase of Bi2W06- 



r- (A u ) r+ (A g ) r 2 - (b u ) r+ (B g ) 



W8i 


87 


54 


63i 


61 


91 


82 


178 


278 


122 


96 


241 


295 


160 


144 


303 


336 


167 


194 


378 


408 


192 


209 


455 


653 


272 


280 


456 


654 


298 


309 






333 


352 






380 


353 






404 


358 






441 


444 






468 


473 






497 


624 






642 


787 






663 


876 






665 








750 








895 







TABLE XII. Frequencies (in cm x ) of the phonon modes at T in the tetragonal IA/mmm phase of Bi2WOe- Values between 
brackets are those previously reported in Ref. |16| . 



r 5 " {E u ) 


r 5 - (E u )[i6] 


It (E g ) 


r 3 " (A 2u ) 


r 3 " (A 2u )[16] 


F+ (Ai 9 ) T 4 - (B 2u ) F+ (B lg ) 


im 


(271£) 


71 


89 


(51) 


176 312 437 


29 


(34*) 


281 


348 


(337) 


767 


137 


(133) 


517 


469 


(437) 




359 


(211) 




674 


(630) 




650 


(525) 
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